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We present a generalized energy-depot model in which the conversion rate of the internal energy 
into motion can be dependent on the position and the velocity of a particle. When the conversion 
' rate is a general function of the velocity, the active particle exhibits diverse patterns of motion 

including a braking mechanism and a stepping motion. The phase trajectories of the motion are 
investigated in a systematic way. With a particular form of the conversion rate dependent on the 
^ l' position and velocity, the particle shows a spontaneous oscillation characterizing a negative stiffness. 

U \ These types of active behaviors are compared with the similar phenomena observed in biology such 

^0 ■ as the stepping motion of molecular motors and the amplification in hearing mechanism. Hence, 

our model can provide a generic understanding of the active motion related to the energy conversion 
and also a new control mechanism for nano-robots. We also investigate the noise effect, especially 
on the stepping motion and observe the random walk-like behavior as expected. 
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^ ' PACS numbers: 87.10.-|-e, 05.45.-a, 47.35.-i 



I. INTRODUCTION 

o 

' ^ , Active motion is a phenomenon found ubiquitously in nature ranging from the colony organization of microorganism, 
J>-^' biological springs, molecular motors, nano-robotics [U, H, H, 0, Q , a flocking of fish and birds to swarming of small 
, insects . All of these types of active motions need processes of energy supply, conversion to motion, and resupply 
for repeated motions. To describe these processes, an energy depot model was introduced by Schweitzer et. al |8|, l9|. 
In the energy depot model, supplied energy is able to induce effectively a negative friction in a certain range of velocity 
yielding an active motion. This model was successfully appHed to a wide variety of active motions [sl. R [sl lol. [Tol [Til. [l2| 
^ mainly assuming that the energy conversion rate into the motion depends only on the quadratic form of the velocity. 
This success of the quadratic energy conversion rate model immediately raises the possibility that other diverse active 
motions can be induced, when the energy conversion rate contains various terms in velocity including the quadratic 

• ' Another motivation for the present study originates from the fields of nano-robotics and artificial molecular mo- 
tions Rapid developments in these area require diverse controlling technology for energy conversion for intended 

purpose. Although it is still at a primitive stage at present, control of the energy conversion will become fine-tuned 
in the near future so that any desired motion could be selectively achieved. Motivated by these considerations, we 
consider various forms of the conversion rate of the internal energy into a motion in this paper. Hence our work 
provides a correlation between the form of energy conversion and the active motion. 

In the following sections, we introduce our generalized energy depot model in detail and consider various forms of 
5J] the conversion rate which could induce diverse active motions including a braking mechanism applicable to nonlinear 
compressive behavior of basilar membrane in the inner ear |,14] . a directed stepping motion shown in the molecular 
motor systems [l^, [l^, [l3| , or a phenomenon of negative stiffness for oscillatory system such as the hair bundle in 
mammals PI- 

Concluding remarks will follow. 
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II. GENERALIZED ENERGY DEPOT MODEL 



A Brownian particle moving in an external force f{x) is governed by the Langevin equation, 

= -fiov + fix) + ^2kBTfioC{t), (1) 

where m is the mass, v the velocity of the particle, /ip the friction coefficient, T the temperature, and the white noise, 
satisfies (C(^)) — and {Ci{t)Cj{t')) — dijd{t~t'). The passive motion can be converted to the active motion when 
the internal energy depot is introduced. The internal energy of the depot, e{t), is regarded as an additional degree of 
freedom for each particle. The energy balance equation for the depot can be described as @ 

a^^q{x)-ce{t)-d{x,v)e{t), (2) 

where a is the timescale of relaxation of the depot, q the rate of the energy influx to the depot, c the rate of energy 
dissipation of the depot, d{x, v) the conversion rate of the internal energy into motion, a = 1 means that the depot 
reacts with a time lag and a that the depot adapts very fast (adiabatic approximation). We note that the 
conversion rate of the internal energy into motion can be, in general, a function of both velocity and position of the 
particle, which contains various combinations of the two variables [a, [l^, [l^l . In order to understand the roles of 
individual terms, we expand d{x, v) as follows, 

oo 

d{x,v) = ai,jx\i^ . (3) 

We can choose ao,o = since the effect of constant conversion rate can be merged into the constant dissipation rate. 
Consequently, the Langevin equation for an active Brownian particle is written as 

dv 



m— = -fJ-oV + fix) + Factive + \/2kBT flQC{t) , 

= -fi{x,v)v + f{x) + V2fcBTK(i), 

where Factive = d{x,v)e{t) /v and ii[x,v) = /io — d{x,v)e{t) /v"^ . The above equations are intuitive in a sense that the 
Brownian particle is governed by an extra active force or that the friction is modified to depend on the space and 
velocity, when external energy is supplied. It is noticeable that the effective friction can be negative or larger than 
the normal friction /xq. The effective friction will be discussed in the next section in detail. 

In this paper, q and c are assumed to be constant for simplicity and main focus is given on the effect of the 
conversion rate. For this purpose, we consider the conversion rate d{v) depending only on the velocity first and c?(x, v) 
in a special form later. Even when only velocity-dependent conversion rate is considered, diverse properties can be 
discussed such as the motion with a braking mechanism, the stability of fixed points, and a stepping motion. A special 
form of d{x, v) is chosen to describe a possibility of a negative stiffness. 



III. MOTION WITH A BRAKING MECHANISM 



The active Brownian particle with d{v) ^ v"^ has been previously treated in detail and it is shown that the friction 
can be effectively negative when the particle moves slowly and increases only up to /xo as the particle moves faster jsl, @1 . 
Hence, for an active particle with d{v) ~ w^, speed can be increased without any limit. Since such an excessive speed 
can damage living organisms, it may be possible that a living organism is equipped with a protective mechanism to 
prevent damage from excessive movement or energy pumping, especially in underdamped oscillating systems. Also, 
in nano-robotics, adoption of this type of protective mechanism may be not only helpful to control the movement, but 
also essential to safeguard the mechanism from overdriving. Thus, study on contributions from higher order velocity 
terms are highly desirable. 

Motivated from these arguments, we first study a symmetric form including up to a fourth order term in the 
conversion rate of the internal energy into motion. Hence, 



d{v) = ao,2V^ + ao^iv'^, 



(5) 
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FIG. 1: (color online) An illustration of the effective friction as a function of the velocity. Solid line shows the effect of the 
fourth order term in the conversion rate, compared with the Schweitzer's result in dashed line, q = 10, = 2.0, Vc = 5.0 and 
0.0.2/ c = 1 are used for the numerical calculation. 



where the positive (negative) 00,4 increases (decreases) the kinetic energy of the particle. Since we are interested in 
the braking mechanism, only the negative case will be treated. By introducing a critical velocity, Vc — \/ao,2/|ao,4|j 
the conversion rate is rewritten as 

d{v) = ao,2v' (1 - ^) • (6) 

It implies that the energy depot reabsorbs the kinetic energy when the velocity goes beyond the critical velocity. This 
mechanism is analogous to the regenerative-brake system in automobiles and electric vehicles [2l|. It is shown that 
this regenerative-braking can not only control the motion more effectively, but also save energy for operation. Thus, 
we name this fourth-order model as a Brownian energy depot model with a braking mechanism. 

To describe the effect of the fourth order term in the conversion rate, let us consider an adiabatic approximation, 
a ^ 0, in which the energy depot adapts very fast. The adiabatic solution of the internal energy of depot yields _£] 

c + ao,2 (1 — w j ^ 
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Hence the effective friction coefficient becomes 



^'^ c + ao,2(l — w^/Wc)w^ 



Note that /i(wc) = A'o- The effective friction describes three different situations ; (i) pumping where [i(v) < 0, (ii) 
dissipation where < iJ,{v) < /zq, and (iii) protection where /xq < /^^(w). The particle cannot speed up when the velocity 
exceeds the critical velocity in this model, since the fourth order term in the conversion rate induces sufficiently large 
friction to the particle. Superficially, the present model appears to be similar to the Rayleigh model [l^l, but with 
different interpretations. 

We believe that the newly introduced braking mechanism may offer a plausible way to discuss any underdamped 
motions and protection from any overreacting motion in living organisms and nano-robotics. Recently by adopting 
this mechanism to the active oscillation of the basilar membrane (BM) in the mammalian ear, we could successfully 
explain the known experimental results and the noise amplification in cochlea p^ . 



IV. STABILITY OF FIXED POINTS 



So far, we have restricted the conversion rate d(v) to be symmetric in the velocity. However, the active particles 
generally possess polarities in their motions [ij, [Tsl . [l6l . [l7| . The polarity of the motion is originated from the 
external force or/and the asymmetric conversion rate. To discuss the generic polarity of the motion, in this section, 
we consider a general form of d{v), 

d{v) = ao^iv + ao^2v'^ + aQ,3V^ + aoAv'^ + ■ ■ ■ . (9) 

It is expected that inclusion of these odd terms allows the active particles perform polarized motions as observed in 
molecular motors such as kinesin and dynein. Now, it is worth to analyze the deterministic dynamics in the {v, e} 
phase space which is governed by 

fJ-o , f , d{v) 

v = V + — + e, 

m m mv (,iUj 

e = q — ce ~ d(v)e, 



where the depot takes time to be filled with energy (a = 1). It is well known that long time behavior of a two- 
dimensional continuous dynamical system possesses only fixed points and limit cycles [1^ . Although it is not possible 
to describe the solution analytically, local stability for existing fixed points can be analyzed through the Jacobian 
matrix and its eigenvalues [l^ . In the following, we discuss the stability of the fixed points obtained numerically 
for various values of a^j. For the numerical calculation, m = 1, q = 10, c = 0.01 and /xq = 20 are used. 

Table U shows the bifurcation with the variation of ao.3 when we consider the conversion rate only up to the third 
order of velocity and there is no external force. There is always at least one stable fixed point or stable focal fixed 
point in the {v, e} space. The trajectories initiated around the stable fixed point or stable focal fixed point are finally 
localized at this fixed point. When —0.16 < ao,3 < 0.16, there are two stable focal points. On the other hand, a limit 
cycle appears if 09,3 > 0.16. In Fig. [21 wc show several phase trajectories on the phase space, {v,e} for several values 
of ao,3 which describe different classifications of fixed points as listed in Table [D Here, we note that overall behavior 
of this system is not much sensitive to variations of ao,i. 

To see the pattern of motion for the limit cycle, we replot the phase diagram for ao,3 = 1.2 in Fig. [3]Ja) with 
corresponding v{t) and x{t) in Fig. [3]Jb) and (c), respectively. It shows that the quasi-triangle shaped hmit cycle 
corresponds to a step-by-step movement of the particle (solid line) while the stable fixed point does to a smooth 
movement (dashed line). The direction of motion depends on the initial velocity. We will discuss this stepping motion 
in detail later because it may be applicable to the molecular motor systems in biology. 

Now, we include the fourth order term in the conversion rate, d{v) and analyze the classification of the fixed points 
as ao,4 is varied. 

d{v) = aoAV + ao,2W^ + 00,3^^^ + ao,4i'''- (11) 

For the numerical analysis, we use ao.i = 0.0002, 00,2 = 2.0, ao,3 — 1.2, and / = 0. We will discuss the effect of the 
external force later. 
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TABLE I: Bifurcation with a variation of ao,3, when d{v) — ao.it' + ao,2«^ + ao,3^;^ ao,i = 0.0002, ao,2 = 2.0, and / = 0. (+/-) 
denotes positive/negative velocity. 

Range of ao.a Classification of the fixed point (s) 

ao,3 < —0.93 • one unstable focal fixed point (-) 

• one saddle point (-) 

• one stable fixed point (+) 

—0.93 < ao,3 < —0.16 • one unstable focal fixed point (-) 
one saddle point (-) 
one stable focal fixed point (+) 

one stable focal fixed point (-) 
one saddle point (-) 
one stable focal fixed point (+) 

one stable focal fixed point (-) 
one saddle point (-) 
one unstable focal fixed point 
with a limit cycle around it (+) 

one stable fixed point (-) 
one saddle point (-) 
one unstable focal fixed point 
with a limit cycle around it (+) 



In Table [III the bifurcations and corresponding nature of fixed points for different values of 09,4 are listed. When 
ao.4 < —0.69, the particle moves smoothly regardless of the initial velocity. As 00,4 increases, the motion starts to 
depend on the initial state of the particle. When —0.69 < 00,4 < 0.92, the particle shows a directed stepping motion 
or a smooth motion depending on the sign of the initial velocity. Interestingly, when ao,4 > 0.92, the particle moves 
stepwise in the positive direction eventually regardless of its initial velocity. We plot this behavior in Fig HI when 
aoA = 1-0. 

So far, we have not included any external forces to the system. However, it has been shown that the active particle 
is able to show a directed motion, without any external forces, depending on the structure of the conversion rate of 
the internal energy into the motion. Another interesting feature is that the directed motion could be stepwise. Such 
a directed stepwise motion appears in the motion of molecular motors in biology. The molecular motors usually carry 
external loads dictated by their functions. In order to simulate this load-carrying capability of the molecular motors, 
we now apply an external force, /, to the active particle. To analyze the bifurcation as / varies, we assume that / 
remains constant during the motion and use that Qq i — 0.0002, ao,2 = 2.0, 09^3 = 1.20, and ag 4 = 1.0. Table Hill 
shows the classification of the fixed points for various ranges of /, where a positive force pulls the particle to the 
positive direction and a negative force acts as a load against the positive direction of motion. 

When the load is large enough, / < —0.89, a limit cycle does not appear, yielding a smooth movement pulled back. 
If —0.89 < / < —0.85, the pulled back particle moves stepwise except when it is initialized close to the left fixed point. 
On the other hand, the particle is pulled back due to the load and moves stepwise regardless of the initial state if 
-0.85 < / < -0.0021. Note that the the saddle point moves from v >0 (when / < -0.2) to w < (when -0.2 < /) 
by crossing the value, v = when / = —0.2. In a very narrow range, —0.0021 < / < —0.0019, the particle is able to 
move stepwise in both negative and positive directions depending on the initial state. When —0.0019 < / < 5.38, the 
particle can overcome the load and move stepwise to the opposite direction of the load. For a sufficiently large force 
to the positive direction, the particle is dragged smoothly to the same direction of the force regardless of the initial 
state of the particle. We plot these trajectories in Fig. O 



-0.16 < ao,3 < 0.16 



0.16 < ao,3 < 0.93 



ao,3 > 0.93 
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FIG. 2: Trajectories in the phase space, {v,e} for several values of ao,3, when ao,i — 0.0002, ao,2 ~ 2.0, and f — 0. Note that 
a limit cycle appears only if ao,3 > 0.16. 



V. STEPPING MOTION 



The limit cycle in our phase space, {v,e} indicates a stepwise motion in time. Moreover a particle carrying a load 
(/, in this model) can perform a directed step-by-step motion overcoming the load. This motion is similar to that of 
a processive molecular motor in a cell. Here, we analyze an engine mechanism of the limit cycle and relate it to the 
walking mechanism of the molecular motors. A noise effect to the stepping motion in our generalized energy depot 
model is also discussed. 

Fig. [5] shows a typical pattern of the limit cycle appearing in our model. This triangle-like limit cycle describes a 
cycling motions between three states which are resting (e > 0, w = 0), accelerating (e < 0, w > 0), and decelerating 
(e ~ 0, w < 0). The energy is supplied during the resting state, consumed for the acceleration, and barely changed 
during the deceleration. Hence, cyclic repetition of this motion leads to the stepping motion as mentioned in the 
previous sections. 

To explain the walking mechanisms of molecular motors such as kinesin and myosin V, a three-state model has 
been proposed 0,12^. The three-state model describes "ATP-binding" to the motor at rest, "ATP-hydrolysis and 
movement" , and "attachment and ADP-releasing" . The "ATP-binding" plays a role of the energy input (e > 0) 
while stalling {ii = 0) and the "ATP-hydrolysis and movement" corresponds to the energy consumption (e < 0) for 
the acceleration {v > 0) in our model. When ADP is released, there is no energy cost. Hence the "attachment and 
ADP-releasing" state corresponds to decelerating state {e = 0,v < 0) in our model. The difference between the 
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FIG. 3: (color online) (a) Trajectories in the phase space, {^ije} when ao,i = 0.0002, ao,2 ~ 2.00, ao,3 = 1.2, and / = 0. (b) 
When the initial velocity is in the positive (negative) direction, the velocity has a steady jerky (smooth) pattern in time, (c) 
Displacement of the particle is step-like (smooth) for the limit cycle (stable fixed point). 
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FIG. 4: (color online) Trajectories when ao,i — 0.0002, ao,2 ~ 2.0, ao,3 = 1.2, ao,4 — 1.0, and / = 0. (a) The phase diagram 
of the particle in the {v, e} space. As shown in Table HH There is one limit cycle around the fixed point with the positive 
velocity, (b) The velocity is jerky regardless of the direction of the initial velocity, (c) The particle shows stepping motions to 
the positive direction even for the initial motion of negative velocity. 
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TABLE II; Bifurcation with a variation of ao,4 when d{v) = ao.it' + ao,2t'^ + ao.3t'^ + cio,4«*, ao,i = 0.0002, ao,2 = 2.0, ao,3 ~ 1-2, 
and f — 0. (+/ — ) denotes positive/negative velocity. 



Range of ao,4 Classification of the fixed point (s) 

aoA < —1.93 • one stable fixed point (-) 

• one saddle point (-) 

• one stable fixed point (+) 

— 1.93 < ao,4 < —0.69 • one stable fixed point (-) 

• one saddle point (-) 

• one stable focal fixed point (+) 

—0.69 < ao,4 < 0.25 • one stable fixed point (-) 

• one saddle point (-) 

• one unstable focal fixed point 
with a limit cycle around it (+) 

0.25 < ao,4 < 0.92 • one stable focal fixed point (-) 

• one saddle point (-) 

• one unstable focal fixed point 
with a limit cycle around it (+) 

0.92 < ao,4 < 2.89 • one unstable focal fixed point (-) 

• one saddle point (-) 

• one unstable focal fixed point 
with a limit cycle around it (+) 

ao,4 > 2.89 • one unstable fixed point (-) 

• one saddle point (-) 

• one unstable focal fixed point 
with a limit cycle around it (+) 



three state model for the molecular motor and our generalized energy depot model is that the molecular motors are 
walking on tracks such as microtubule or actin filaments. However, the interaction between the motor protein and 
the microtubule or the actin filament can be considered in our model by introducing the corresponding external force. 
We will report on this work later elsewhere. 

So far, we have neglected the stochastic noise contribution to the motion of the active particles. However, the 
thermal fluctuation is inevitable and expected to play an important role especially in the motion of the molecular 
motors and nano-robots. Hence we include the stochastic term with the noise strength \/2fcsT'/^o in calculation 
of the motion. We calculate the effect of the stochastic noise to the stepping motions described in Fig. Eljd) and 
(e) in which the applied loads are / = —0.002 and / = —0.001, respectively. The stochastic stepping motions are 
plotted in Fig. [71 when the strength of the noise is 0.6. As expected, the particle can move forward or backward 
as a consequence of the noise fluctuation, yielding irregularly directed stepwise motion as observed in the motion 
of the molecular motors llSl Il6l. llTI . Here, we compare the present work to the existing studies on the directed 
motion of active particles [id . I 111 . Il2 |. In those works, the asymmetry of space was described by an external rachet 
potential [l^, [O, However, in our model, this asymmetry is given by the conversion function, i.e., the en ergy 
depot itself. In our model, the noise introduces occasional backward motion as observed in experiments p^ . [igI . I17| |. 
unlike previous models. Therefore, further study on the effect of the external ratchet potential on the present model 
appears necessary to understand the physical origin of the motional polarity. 
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TABLE III: Bifurcation with a the variation of the external force when ao,i = 0.0002, ao,2 = 2.0, 00,3 = 1.2 and ao,4 = 1.0. 
(+/— ) denotes positive/negative velocity. 



Range of / 


Classification of the fixed point (s) 


/ < -2.77 


• one stable focal fixed point (-) 

• one saddle point (+) 

• one unstable fixed point (+) 


-2.77 < / < -0.89 


• one stable focal fixed point (-) 

• one saddle point (+) 

• one unstable focal fixed point (+) 


-0.89 < / < -0.85 


• one stable focal fixed point 
with a limit cycle around it (-) 

• one saddle point (+) 

• one unstable focal fixed point (+) 


-0.85 < / < -0.2 


• one unstable focal fixed point 
with a limit cycle around it (-) 

• one saddle point (+) 

• one unstable focal fixed point (+) 


/ = -0.2 


• one unstable focal fixed point 
with a limit cycle around it (-) 

• one saddle point at the resting state {v = 0) 

• one unstable focal fixed point (+) 


-0.2 < / < -0.0021 


• one unstable focal fixed point 
with a limit cycle around it (-) 

• one saddle point (-) 

• one unstable focal fixed point (+) 


-0.0021 < / < -0.0019 


• one unstable focal fixed point 
with limit a cycle around it (-) 

• one saddle point (-) 

• one unstable focal fixed point 
with a limit cycle around it (+) 


-0.0019 < / < 5.38 


• one unstable focal fixed point 

• one saddle point (-) 

• one unstable focal fixed point 
with a limit cycle around it (+) 


5.38 < / < 6.496 


• one unstable focal fixed point (-) 



• one saddle point (-) 

• one stable focal fixed point (+) 



/ > 6.496 



one unstable fixed point (-) 

one saddle point (-) 

one stable focal fixed point (+) 
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FIG. 5: (color online) Trajectories in the phase space, {v, e} for several values of / chosen from Table Hill when ao.i = 0.0002, 
ao,2 = 2.0, ao,3 = 1.2, and ao,4 = 1.0. (a) The particle is pulled back in a smooth way due to a heavy load.(b) The pulled 
back particle moves stepwise except when it is initialized close to the left fixed point, (c) The particle is pulled back due to the 
load and moves stepwise regardless of the initial state, (d) The particle is able to move stepwise in both negative and positive 
directions depending on the initial state, (e) The particle can overcome the load and move stepwise to the positive direction, 
(f) A sufficiently large force to the positive direction drags the particle smoothly to the positive direction regardless of the 
initial state of the particle. 



VI. NEGATIVE STIFFNESS 

The conversion rate of the internal energy into the motion is in general a function of space and velocity [l^l • As 
we have observed that the velocity dependent conversion rate can induce an active behavior in certain velocity range, 
it is expected that a position dependent conversion rate may also introduce an active behavior in oscillation. In this 
section, we consider a special form of the space-dependent conversion rate to show that our generalized energy depot 
model induces an active behavior in certain spatial range. This active behavior in an oscillatory system appears as 
the negative stiffness of membrane or bundle. As an example, it is known that the spontaneous oscillation of the 
mechanosensitive hair bundle in the inner ear is strongly related to the negative stiffness of the bundle [Ts] . Therefore, 
we believe that it is desirable to investigate contributions from a conversion rate function, which depends on velocity 
and position simultaneously. 
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FIG. 6: The engine mechanism of the limit cycle. The motion can be approximately divided into three states; resting (e > 0, 
V = 0), accelerating (e < 0, ?) > 0), and decelerating(e ~ 0, ?) < 0). 



Motivated by this observation, we consider the following form of the conversion rate, 

d{x,v) — ai^ixv + as^ix^v. (12) 
For the stationary state, e = 0, the equation of motion becomes 

mv = -fiov - Kc«x + /, (13) 

where the effective stiffness is 

KcS = K - {ai,i + a3,ix'^)e. (14) 

To avoid infinite amplitude of vibration, 03^1 should be negative. Note that if the second term exceeds the physical 
stiffness k, KeS could be negative in the region (ai,ie — «;)/|a3,i |e < x < ^ (ai^ie — K)/|a3,i|e. Fig. [8] plots the 
extention force KoffX versus x. This type of the negative stiffness has been experimentally observed and explained 
using a two-state model [3]. Thus, it will be interesting to investigate possible connections between the position 
dependent conversion rate and the two-energy-state model. Another possible application of the position dependent 
conversion rate is introduction of a protective mechanism for the system from excessive displacement. We believe that 
a suitable choice of the conversion rate can provide such a protective mechanism for oscillatory nano-robots. 
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FIG. 7; (color online) The stochastic stepping motion of the active particle with a constant load / when ao,i — 0.0002, 
ao,2 = 2.0, ao,3 = 1.2, ao,4 = 0.93, and the noise strength is 0.6; (a) / = -0.002 (b) / = -0.001. 



VII. CONCLUSION 



In this paper, we have presented a general model of active Brownian particle in which the conversion rate of the 
internal energy into the motion is a general function of the velocity of the particle. We have shown that when the 
conversion rate depends only on v"^ and v^, an active amplification and a protective braking mechanism appear. We 
believe that this active amplification with a braking mechanism can provide a new paradigm for any underdamped 
biological systems and nano-machines. It is conceivable that the evolutionary tactic in biological system may have 
provided such protective mechanisms to any underdamped biological units to protect them from excessive oscillations. 
Such a self-adapting mechanism will be also needed in designing nano-machines. 

When the conversion rate is asymmetric in velocity, various interesting behaviors appear. One of the most re- 
markable behaviors is a directed motion without any external forces. Such a directed motion suggest a possibility of 
overcoming any extra load on the particle. Another most remarkable behavior is a stepping motion described as a 
limit cycle. Such a stepping motion has a polarity hence being applicable to the molecular motors. 

When the conversion rate depends on not only the velocity but also the position of the particle, it has been shown, 
for a very special case, that the effective stiffness of an oscillatory system could be negative. Such a behavior has 
been observed in the hair bundle motion in the inner ear. Hence it also suggests a possible connection between our 
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FIG. 8: Extension force Hesx versus x when d{x,v) — ai^ixv + as^ix^v. In this plot, we use that ai,i = 2.1, a3,i — —1, k = 1, 
and e — 1. 

generalized energy-depot model and a relevant biological system. 

We have shown that various interesting active motions which exist in real world can be induced for active Brownian 
particles by considering various forms of the energy conversion into motion. It is quite remarkable that our general 
energy-depot model could describe such a variety of motions even though the generic relations between our model 
and the real world are to be investigated further. 
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